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ABSTRACT 

■ In light of recent Super-Kamiokande data and global fits that seem to exclude both pure 

G^ ■ i^fi —>■ i^s oscillations of atmospheric neutrinos and pure Ue — > Ug oscillations of solar neutrinos 
^ ! (where Ug is a sterile neutrino), we reconsider four-neutrino models to explain the LSND, 

OO I atmospheric, and solar neutrino oscillation indications. We argue that the solar data, with 



O 



X 



the exception of the CI results, are suggestive of — > Us oscillations that average to a 



O ' probability of approximately |. In this interpretation, with two pairs of nearly degenerate 
mass eigenstates separated by order 1 eV, the day-night asymmetry, seasonal dependence, 
Ph, and energy dependence for ^B neutrinos should be small. Alternatively, we find that four- 
Qh! neutrino models with one mass eigenstate widely separated from the others (and with small 
^ I sterile mixings to active neutrinos) may now be acceptable in light of recently updated LSND 
results; the ^''Cl data can be accommodated in this model. For each scenario, we present 
simple four-neutrino mixing matrices that fit the stated criterion and discuss future tests. 



1 Introduction 

Accelerator, atmospheric, and solar neutrino data give evidence for neutrino oscillations 
and thus for neutrino masses and mixing. The LSND accelerator experiment finds a small 
~^ appearance probability and a mass-squared difference Sml^^^) > 0.2 eV^ [|T], |^. The 
atmospheric experiments measure the —>■ survival probability versus both path-length 
L and neutrino energy E,,. The amplitude for Vx {x ^ e) oscillations is inferred to be 

maximal or near-maximal with (5m^rpj^ ~ 3 x 10~^eV^ p, ^ The combined solar neutrino 
experiments determine the z/g — > z/g survival probability to be ~ 0.3 to 0.7, depending on 
the neutrino energy 0, |^, ^ and ^w^solar ~ lO '^ required for consistency with the 



null result for the Ve ~^ probability measurement in the CHOOZ experiment [|T0[. Thus, 
taken together, the data require three distinct 5m^. With three neutrinos, there are only 
two independent Sm'^, so a fourth sterile neutrino (us) needs to be invoked in addition to Ue, 
Vai and Vt- Having no weak interactions, the sterile neutrino escapes the A,^ ~ 3 constraint 



from the invisible decay width of the Z-boson 



Successful global descriptions of the oscillation data have been made in a four-neutrino 



framework |T^, |T^, |l^, |T6[. Compatibility of the LSND result with null results of acceler- 



ator and reactor [|I^ disappearance experiments was found |T2| , to favor a 2 -|- 2 mass 
spectrum, with the mass-squared differences 5w^solar ^^atm separated by the larger 
difference 5?7t.lsnd5 over a H-3 spectrum, with one mass eigenstate widely separated from the 
others. In the simplest of the 2 + 2 models, the oscillations are Vj. (atmospheric) and 

(solar), or alternatively Vg (atmospheric) and (solar). More generally, 

sterile neutrinos may be involved in both atmospheric and solar oscillations. In the 2 + 2 
mixing schemes, the sterile flavor content must be significant in the solar or atmospheric 
oscillations or in both. 

New results from Super-Kamiokande (SuperK) 0, 0] impact critically on oscillations to 
sterile neutrinos. The zenith angle dependence of high-energy atmospheric events, along with 
the neutral current vr" production rate, exclude pure oscillations at 99% C.L. The 

electron energy dependence of solar neutrino events, along with the absence of a significant 
day-night effect, exclude pure z/g oscillations at 99% C.L. when all of the solar data, 

including the ^^Cl results, are taken into account. Thus a dominant involvement of Vs in 
either atmospheric or solar oscillations is brought into question. Also, LSND has reported 
new results |^ with a somewhat lower average oscillation probability than before and the 
KARMEN experiment excludes part of the LSND allowed region. Consequently, previous 
exclusions [ll2| , |l^ of the 1 + 3 schemes need to be reassessed. In the context of this new 
data, we reconsider four-neutrino oscillation models and the implications for the existence 
of a sterile neutrino. 



2 2 + 2 Models 

The 2 + 2 models have two nearly degenerate pairs of neutrino mass eigenstates giving 
(5m^rpj^ and 5^«solar separated by the LSND scale 5?71lsnd ~ 1 ^V^. Maximal z/3, z/2 mixing 
describes the atmospheric data and maximal vq, vi mixing describes the solar data. The mass 
hierarchy may be normal {m^ ^ mi) or inverted {m^ -C mi). The ^^t^solar splitting must 
be smaller than 10~^ eV^ to avoid the CHOOZ reactor constraint on Uf. — >■ Uf, oscillations 
and larger than the ~ 10~^° eV^ of just-so vacuum oscillations PO] so that the oscillations 



average to give an approximately energy-independent solar Ve survival probability. 
2.1 Solar Neutrinos: maximal Uf, Ug oscillations 

The new SuperK solar neutrino data show a remarkably fiat recoil electron energy distribu- 
tion from 5 MeV to 14 MeV, with average value 

data/SSM = 0.465+°:°}^ , (1) 

where SSM is the Standard Solar Model prediction [Q. The flux normalization in the SSM 
is somewhat uncertain, and the above result is suggestive of rapid oscillations with maximal 
amplitude that give an average oscillation probability, {P{i^e ^e)), of approximately |. 

Regeneration of in the Earth, due to the effects of coherent forward scattering of 
the neutrinos with matter pl[, causes a day-night variation in the measured neutrino flux. 
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The smallness of the day-night effect 2{D - N)/{D + N) = -0.034 ± 0.022 ± 0.013 [0] in 
the SuperK data (1.3(7 from zero) puts strong constraints on solar solutions with matter 
effects [0, |2^, It has been noted that there can be a day-night effect even for 



maximal mixing |]23|, but calculations for active neutrinos show that it is small (~ 1- 
2 %) for ^m^Qia^j, > 5 X 10"^ eV^ or 5ml^^^-^ < 10^"^ eV^ 0. In z/g z/g oscillations, the 
difference in coherent scattering amplitudes in matter is \/2Gp{Nf. — A^„/2) \f2GpN(,l2 
(where and iV„ are the electron and neutron number densities, respectively, and N^. ~ 
in the Earth). The corresponding amplitude difference for — (or —> u^-) is \/2GFNe. 
Therefore, for maximal mixing the day-night effect should be smaller with a sterile neutrino 
than with an active neutrino, but it would not be zero. 

The regeneration effect in the earth for neutrinos detected at night raises the value of 



{P{ye i^e)) for maximal mixing with active neutrinos from 0.50 to about 0.54 ||2^. For 
sterile neutrinos, we roughly estimate that {P{ve —>■ ^e)) — 0.52 in SuperK for maximal 
mixing with regeneration. Regeneration also causes a weak dependence on energy of the 
solar neutrino suppression. For maximal mixing of v^. with an active neutrino, {P^i^e ^e)) 
is slightly higher at higher neutrino energies consistent with the weak trend of the 
SuperK data; less energy dependence should be present with z/g oscillations to Vg- 

Assuming (-P(z/e — ^ Ve)) = 0.52 for maximal i/g oscillations, the SuperK result in 

Eq. (|l]) would be reproduced by a *B flux normalization of n = 0.89. The SSM predictions (in 
SNU) for the Gallium experiments are (70, pp), (34, ^Be), (3, pep), (10, CNO), and (12, ^B), 
giving a total of 129 SNU; the observed values are 75.41^;^ (SAGE), 77.5l^| (GALLEX) and 
65.8l|oi (GNO). With the above *B normalization and (-P(z^e ^e)) = |, the predicted rate 
in the Gallium experiments is 63.8 SNU. Thus all of the Gallium measurements are consistent 
within 2(7 of the value that would be found for maximal amplitude — i^s oscillations that 
average to \. The new data from GNO are especially suggestive of this interpretation. A 
least-squares fit to the SuperK and Gallium rates with the ®B flux normalization n as a free 
parameter gives a best fit value oi n = 0.90 with = 5.5 for 3 degrees of freedom, which 
corresponds to a 14% goodness of fit. In this calculation we have not taken into account any 
regeneration effect for the Gallium data, which may improve the fit. 

The discrepant measurement in the above interpretation is the ^"^Cl value of data/SSM = 
0.33 ± 0.03 from the Homestake mine experiment, which would require a significant energy 
dependence of the solar flux suppression. Most global fits to neutrino oscillation data 
include the ^''Cl data (and often disregard the LSND data) and then z/g oscillations are 
excluded. We instead suggest the possibility that the solar data/SSM flux ratio is relatively 
flat over the entire 0.233 MeV to 14 MeV energy range and is described by Ve i^s oscillations 
with maximal mixing. 

An important test of approximate constant suppression of the solar neutrino spectrum 
will be the BOREXINO [^] experiment, which can measure the ^Be component of the solar 
neutrino flux. Nearly complete suppression of the ''Be component is needed for the ^^Cl data 
to be consistent with the suppression of the ^B component measured by SuperK. However, 
a suppression of the ^Be flux similar to the data/SSM measured by SuperK would favor 
maximal Ve mixing of solar neutrinos. The SNO p6|] and ICARUS |^ experiments 

will also provide critical tests of the *B flux suppression, and SNO will test maximal z/g 
Vs oscillations through the NC/CC ratio (which should be the same as the value without 
oscillations) . 
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2.2 Mixing matrix for the 2 + 2 model 



We assume a 2 + 2 scenario in which one pair of nearly degenerate mass eigenstates has 
maximal z/g —>■ Vs mixing for solar neutrinos and the other pair has maximal (or nearly 
maximal) z/^ — > Vr oscillations for atmospheric neutrinos. Small off-diagonal mixings between 
i^e, T^s and z/^, Vj. can accommodate the LSND results. Neglecting CP-violating phases for 
the present, an approximate mixing with these properties is 
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where l^w^iol ^ I'^'^Ll ^ I ^"^21 1- notation of Ref. 

angles 6101 = 6*23 = vr/4, 6*02 = 6^03 = 0, and dx2 = On 
probabilities are given in Table 0. The parameter e is determined from the LSND data to 
be approximately 



we have chosen the mixing 
with e = sin6'i3. The oscillation 



0.016 eV' 

. I^^lsndI 



0.91 



(3) 



where 5mLSND is restricted to the range 0.2 to 1.7 eV^ by the BUGEY and KARMEN 
experiments; for 5m\^ 



LSND 



6 eV 



0.022 is marginally possible. 

Refs. [g, 0, |ll |1|, 



are also ac- 



Other forms of the mixing matrix (see, e.g. 
ceptable, provided that the Ve-Vs mixing is nearly maximal. Future short and long-baseline 
experiments will be useful in determining the mixing matrix. For example, there are no 
z/g v-r oscillations at short baselines for the mixing in Eq. (^, but there could be if the 
\Je2 and [/e3 mixing matrix elements are not equal in magnitude or have different phases 
Refs. 



[see, e.g. 



2c 



0[), in which case unitarity requires t/ro 7^ and/or \Jr\ 7^ 0. 
Also, both Vf, z/^ and v^, — > oscillations are possible at long baselines, and there can 
be CP violation in these channels if the four-neutrino mixing matrix is not real (see, e.g.. 



Refs. 28 



31, 



3 1 + 3 Models 

Another way to circumvent the SuperK exclusion of sterile neutrinos is to assume that the 
solar and atmospheric oscillations are approximately described by oscillations of three active 
neutrinos; then the LSND result can be explained by a coupling of z/g and through small 
mixings with a sterile neutrino that has a mass eigenvalue widely separated from the others. 
This 1 + 3 model was previously disfavored by incompatibility of the LSND result with null 
results of the CDHS and BUGEY experiments |jl2|, |T^. However, in newly updated results 
from the LSND experiment [Q], the LSND allowed region is slightly shifted, and this opens a 
small window for the 1+3 models. The solar neutrino data can be explained by oscillations of 
three active neutrinos in the usual way; the results of the '^''Cl experiment can be reconciled 
with those of SuperK in part because oscillations of z/g to active neutrinos in SuperK also 
show up as neutral current events (with a rate of about 1/6 of the charge-current rate), 
thereby giving a higher data/SSM in SuperK than in the ^^Cl experiment. 
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By way of example, consider the approximate mixing 
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(4) 



where e and 6 are small and both the flavor and mass eigenstates are reordered to reflect the 
fact that z/q (which is predominantly Ug) is the heaviest state. In the notation of Ref. [28|, 
we have chosen the mixing angles Oqi = 612 = 7r/4, 6*02 = O23 = 0, sin 6^03 = e, and sin^^ia = 6. 
Here the 3x3 submatrix that describes the mixing of the three active neutrinos has the 
bimaximal form [Q. The oscillation probabilities are given in Table For this mixing 
matrix the leading oscillation amplitudes at the ^^lsnd scale in the 1 + 3 scheme are (when 



e,5 < 1) 



4e' 



^BUGEY 

for disappearance in the BUGEY experiment, 

^CDHS ^ 4(5^ 

for Ufj^ disappearance in the CDHS experiment, and 

^LSND = 4e^5^ 



(5) 



(6) 



(7) 



for z/^ 



Vf, appearance in the LSND experiment. Including the subleading oscillation at the 
scale, and assuming that the leading oscillation averages. 



for — >■ disappearance in atmospheric neutrino experiments, where 5m1^^ is in eV^, L 
in km, and E in GeV. The relative normalization of the atmospheric and v^. fluxes is well 



known ; the predicted effective change in this normalization at the 5rn\^^ oscillation scale 



is (1 — 25^) /(I — 2e^). In the region where the CDHS constraint on 5 is weak or nonexistent 
('^"^LSND ~ 0-4 eV^), BUGEY constrains e to be less than about 0.1; conservatively assuming 
that the atmospheric data constrains the relative normalization to within 10%, we find 
5 < 0.24. 

The upper limits on e (from Abugey) and 5 (from Aatm and Acdhs) and the allowed 
values of Alsnd vary with Sm^^^-Q- With the previously reported LSND data |jl|, there was 
no value of ^^tj-lsnd which the constraints on e and 5 were consistent with the LSND 
99% C.L. allowed range of Alsnd; the maximum allowed values of e and 5 always implied 
an upper limit on Alsnd that was below the LSND measured value. However, with the 
recent shift of Alsnd to lower values, there are now three small ^^lsnd islands where the 
LSND 99% C.L. region is compatible with the BUGEY, CDHS, and KARMEN constraints: 



5m 



LSND 



0.9 eV^ and ^^lsnd 



1.7 eV^, for which the BUGEY constraint is somewhat 



less restrictive, and ^^lsnd — 6 eV^, for which both the BUGEY and KARMEN constraints 
are somewhat less restrictive. The E776 experiment at BNL gives a somewhat tighter 
constraint than KARMEN at (5mLSND = 6 eV^, but a small region is still allowed here. 
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Except for a small region near ^^t^lsnd — 0-2 eV^ (see below), the combined data are still 
inconsistent with oscillations having ^?t^lsnd < 0-9 ^^^5 ^^lsnd when the LSND 

90% C.L. allowed region is used. The above analysis is summarized in Table ^; acceptable 
values of ^ttilsnd those for which (4e^5^)niax lies within the range of Alsnd allowed by 
LSND and KARMEN. 

There is no CDHS constraint for ^^lsnd < 0-25 eV^, which suggests that ^^lsnd — 
0.2 eV^ may also be allowed. However, at these 5w^lsnd values 6 must be more than 0.5 to 
reconcile the BUGEY limit with the LSND measurement, which is not consistent with the 
relative normalization of the and Uf, fluxes (see above). 

The model of Eq. (^) has bimaximal mixing in the 3x3 active sector. Similar results for 
the 1 + 3 model can be obtained for any 3x3 submatrix that can describe the atmospheric 
and solar data for active neutrinos, as long as e and 6 are small. This is easily seen by 
realizing that the general expressions for the short-baseline amplitudes in the 1 + 3 model in 
Eqs. (|^)-(0) depend only on the magnitudes of the mixing matrix elements f/go and U^o [e and 
S, respectively, in Eq. (H)]. If Uro 7^ 0, then there can be z/^ — > Ur and z/g — > oscillations 
at the 5w^LSND scale. If Ues 7^ 0, then there can be —>■ and z/g Vt oscillations at the 
5ml^^ scale. 

Because the allowed windows are already severely constrained by accelerator and reactor 
data, future experiments at short baselines could easily test or rule out the 1 + 3 scenarios. 
MiniBooNE will search for oscillations and test Alsnd- Reactor experiments 



such as at Palo Verde [0, KamLAND [Q, and the proposed ORLaND ||3^, will test 
^BUGEY, but only ORLaND is expected to significantly improve the bound on e. A future 
precision measurement of z/^ disappearance at short baselines could test Acdhs- 

4 Discussion and Conclusions 

The recent Super-Kamiokande data and global fits present new constraints on the mixing 
of light sterile neutrinos. It is important to assess the viability of a light sterile neutrino, 
because its existence is required if the solar, atmospheric, reactor, and accelerator data are 
all to be understood in terms of neutrino oscillations. 

We have examined at a qualitative level two schemes for four-neutrino mass and mixing 
which nearly accommodate all present data. In the 2 + 2 scheme, two neutrino pairs are 
separated by the LSND mass scale. One pair maximally mixes z/^ and z/^ to explain the 
atmospheric data, while the other pair maximally mixes and Vs to explain the solar deficit 
with energy-independent oscillations. All data except the low Ve capture rate on ^^Cl are 
explained in this model. A very small day-night effect is expected, at the level of a per cent 
or less. Suppression of the ^Be component of the solar neutrino flux by approximately 50% 
is expected for the BOREXINO experiment. 

In the 1 + 3 scheme, the three active neutrinos are separated by the LSND scale from 
a sterile neutrino state. The active neutrino 3x3 submatrix explains the atmospheric and 
solar data, including the ^''Cl rate. In our example we assigned the bimaximal model [33| to 
this submatrix, but any 3x3 active submatrix that describes the atmospheric and solar data 
may be used. A small mixing of z/g and via the sterile state explains the LSND data. The 
smaller oscillation amplitude recently reported by the LSND collaboration allows marginal 
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accommodation with the null and z/^ disappearance results previously obtained in reactor 
and accelerator experiments. 

The above two classes of models allow a viable sterile neutrino having either large mixing 
with Vf. or small mixing with all active neutrinos. In both classes we made the simplifying 
assumption that Vr-Vf^ mixing is negligible. Present data may allow considerable v^-Vs mixing, 
in which case both solar and atmospheric neutrino oscillations could have a sizable sterile 
component. Quantitative fits to atmospheric and solar data are needed to determine how 
much Vt-Vs mixing is allowed (see, e..g, Ref. |40[|). 

Neutrinos also provide a hot dark matter component, which is relevant to large-scale 
structure formation From the upper limit of = 2.5 to 3 eV on the effective mass 



obtained from tritium beta decay endpoint measurements |^2[| , the mass mmax of the heaviest 
neutrino is bounded by [H3i 



-^"^LSND ~ "^max < \/ m| + ^m^gND ■ (9) 

The contribution of the neutrinos to the mass density of the universe is given by Q^, = 
J2 ^u/ (/i^93 eV), where h is the present Hubble expansion parameter in units of 100 km/ s/Mpc 
@]; with h = 0.65, both the 2 + 2 and 1 + 3 models give fi^ > 0.02. For an inverted neutrino 
mass spectrum (where the Ue is associated predominantly with the heavier neutrinos), the 



bound is y^^lsnd ~ ""^max ~ ^t^, which yields > 0.02 and Q^, > 0.07 for the 2 + 2 and 
1 + 3 models, respectively. The MAP||45[| and PLANCK |^ satellite measurements of the 



cosmic microwave background radiation may be sensitive to these neutrino densities. 

The existence of the LSND mass gap will be tested by the MiniBooNE experiment [36|. A 
definitive test of the presence of the sterile state in the 2 + 2 model will be the measurement 
of the NC/CC ratio for solar neutrinos by the SNO experiment. With the solar solution 
given by maximal u^. — > Us mixing in this model, the NC should show the same suppression 
as the CC. The existence of the sterile neutrino in the 1 + 3 model will not be tested by SNO 
NC/CC data since the ratio is approximately that of I'e oscillations to active neutrinos. The 
sterile state in the 1 + 3 model can be tested by searches for small amplitude oscillations at 
short baselines. 

Observation of the flavor ratio of extragalactic neutrinos may serve as a test of the models. 
Oscillations of neutrinos from distant sources will have averaged, leaving definite predictions 
for flavor ratios. If cosmic neutrinos are mainly produced in pion/muon decay, their initial 
flavor ratio is z/^ : z/^ : z/g : i^s ~ : 2 : 1 : 0. A simple calculation |^ then gives the 



asymptotic ratios 1:1: 0.5 : 0.5 for the 2 + 2 model with two maximally-mixed pairs, and 
1:1:1:0 for the 1 + 3 model with bimaximal mixing of the active neutrinos. 

A virtue of active-sterile neutrino oscillations is that they may aid r-process nucleosyn- 
thesis of heavy elements in neutrino-driven supernovae ejecta. The basic requirement is that 
the Uf. flux be diminished in the region where inverse /5-decay would otherwise transform neu- 
trons into protons. The 2 + 2 and 1 + 3 mass spectra and mixing matrices presented here are 
of the forms previously discussed to enhance r-process nucleosynthesis. In the 2 + 2 model, 
this is accomplished with a two-step process: first, a matter-enhanced v^/vt ~^ i^s transition 
beyond the neutrino sphere removes the energetic v^/fr before they can convert to ener- 
getic z/g, and then a large z/g — >■ z/^/z/^ transition reduces the z/g abundance the required 
mass-squared parameter is 5m\^^-^j) > 1 eV^. In the 1 + 3 scheme, the relevant oscillation 
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is a matter-enhanced — > z/^, obtained with Sm^^^-Q ~ 2 eV^ and sin^ 26*23 ~ 10 ^9[. It 



appears that the r-process enhancements discussed in Refs. P8| , |49|] may be accomphshed by 
tuning the small parameters in our mixing schemes. 
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Table 1: Oscillation probabilities in the 2 + 2 model defined by Eq. (H), to leading order 
for each oscillation scale. The oscillation arguments are defined by A^- = 1.27Sm'^L/ E, with 
in eV^, L in km, and E in GeV. 
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Table 2: Oscillation probabilities in the 1 + 3 model defined by Eq. to leading order 
for each oscillation scale. The oscillation arguments are defined by A^ = 1.27(5m^L/i?, with 
Sm'j in eV^, L in km, and E in GeV. 





e 




T 


s 


e 


1 - 4e^ sin' Alsnd 
- sin^ Asolar 


sin' Alsnd 
+ |sin2 Asolar 


i sin^ Asolar 


4e' sin' Alsnd 
+ 1(^2 _ 2^2) Asolar 


/i 




1 - A5'' sin' Alsnd 
- sin^ Aatm 
- i sin^ Asolar 


sin^ Aatm 
-isin^ Asolar 


4:5' sin' Alsnd 
-5^ sin^ Aatm 
+ I(2e2- 52) sin2 Asolar 


r 






1 - sin^ Aatm 
- i sin^ Asolar 


52 sin2 Aatm 
+ i(2e2- 52) sin2 Asolar 


s 








1 - 4{e' + 5') sin' Alsnd 
-6\2e^ + 6^)sin^AATM 



11 



Table 3: Summary of constraints on four-neutrino oscillation parameters in the 1+3 scheme. 
The upper limit on Alsnd at (J'tilsnd — 6 eV^ is from the E776 experiment at BNL All 
experimental limits are at 90% C.L., except for LSND, which is at 99% CL. 
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